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DATA REDUCTION ANALYSIS AND APPLICATION TECHNIQUE DEVELOPMENT 
FOR AT^WSPHERIC TRACE GAS CONSTITUENTS DERIVED FROM REMOTE 
SENSORS ON SATELLITE OR AIRBORNE PLATFORMS 

By 

Joseph C. Casas^ and Shirley A. Campbell^ 

INTRODUCTION 

This report sunsmarizes work performed under NASA Grant No. NSG-ODS. 
Research efforts concentrated on four comprehensive program areas; 

(1) development of a STAR-100 version of a radiative transfer 
model similar to the SMART program, 

C2) support of the St. Petersburg flight test of MAPS brassboard 
instrument, 

(3) support for the surdner MONEX program, and 

C4) investigative studies for the implementation of the MAPS 
experiment to the OSTA-1 payload for Space Shuttle. 

The primary objective of this project was to investigate the appli- 
cability of the gas filter correlation radiometer (GFCR) to the measurement 
of tropospheric carbon monoxide gas from airborne and earth-orbiting 
platforms. To this end research has foci^sed on the operational application 
of the GFCR technique to the remote measurement of CO from low-altitude 
aircraft platforms. 

An assessment of the GFCR measurement system to a regional measurement 
program was conducted through extensive aircraft flight testing of several 
versions of the GFCR. This test program was conducted by personnel at 
NASA/ Langley Research Center (LaRC) and complemented by an Old Dominion 
University (ODU) comprehensive research program. This research program' 

Research Assistant Professor, Department of Geophysical Sciences, Old 
Dominion University, Norfolk, Virginia 23508. 

^Research Associate, Old Dominion University Reseaich Foundation, P.O. Box 
6369, Norfolk, Virginia 23508. 



involved investigative work in the areas of flight -test planning and 
coordination, acquisition of verifying CO measurements, determination 
and acquisition of supporting meteorological data requirements, and 
development of supporting computational software. 

NOMENCLATURD 

concentration of absorbing gas, ppm 

C8 Boltzmann distribution constant, cm Kelvin 

E monochromatic upwelling radiance, W cm’'‘’sr“^ 

E* energy of the lower sts*,e, cm"^ 

f Chapman function, dimensionless 

h sensor altitude index, dimensionless 

h' uppermost layer altitude index, dimensionless 

H wavenumber-dependent sun irradiance at the top of the atmosphere, 
® W cm’lsr'^cm"^ 

k atmospheric layer .index, dimensionless 

1 thickness of layer, cm 

m constant, 1.5 for water vapor and 1.0 for all other molecules, 

dimensionless 

N° Planck blackbody radia.tion, W cm'^sr”^ 

p mean pressure of layer, atm (1 atm » 1.01325 E+5 N/m^) 
p equivalent pressure, atm 

adjusted spectral line intensity, atm’^cm'^ 

So spectral line intensity, atm'^cm’^ 

T layer temperature, Kelvin 

Tq reference temperature corresponding to spectral line parameters, 
Kelvin 

Tg surface temperature, Kelvin 

z altitude index, dimensionless 



adjusted spectral lino halfwidth, cra"^ 

«0 lino halfwidth, cm”^ 

Lorents line shape, dimensionless 
E wavenumber -dependent surface emittance, dimensionless 
wavenumber-dependent absorption coefficient, atm'^cMi”' 

T gasec^us transmittance at a particular altitude, dimensionless 

0 sun zenith angle, degrees 

0 ) wavenumber (inverse wavelength), cm“^ 

w averaged wavenumber, cm~^ 

Subscripts: 

1 gas species number 

n spectral line number 


STAR- 100 
Introduction 

The use in recent years of high-sensitivity and high effective spectral 
resolution (less than 0.1 cm"^3 instrvimentation, such as MAPS, in the remote 
measurement of trace atmospheric gases has required refinements in the 
computational methods used in calculating the atmospheric transmittance in 
the infrared spectral region. The principal techniques of detecting trace 
atmospheric gases, specifically gaseous pollutants, have been described by 
Ludwig (ref. 1) . The success of these techniques not only depends upon 
the ability of the instrument technique to discriminate between the 
pollutant and interfering spectral lines, but also requires an accurate 
knowledge of high-resolution pollutant and atmospheric spectra, This can 
be accomplished by applying a line-by-line atmospheric radiative transfer 
program. 

The program presented in this report is an improved version of the 
line-by-line program written by Casas and Campbell and described in 
reference 2. The program was written for a vector processing machine, the 
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STAR-100, in a langtiaga (SL/1) specifically designed to take advantage of 
its parallel processor capabilities, The objective* was to develop an 
efficient, generalized, Une-by-line, atmospheric radiative transfer 
program, which would be readily adapted to many different Infrared sensor 
research needs. This report describes the program, the analytical 
concepts upon which it is based, as well as the STAR-100 computer and the 
SL/1 programming language. 

The STAR-100 (ref. 3) is a large-scale, high-speed digital computer 
employing parallel processing and a virtual memory. These features 
permit fast matrix-oriented calculations with practically unlimited 
storage capacity. The STAR-100 is a vector machine; i.e,, the optimum 
performance of the computer occurs with contiguous floating point operands 
defined by a location and length, It achieves a high performance level 
through use of a parallel processor that causes the elements of vectors 
to be fetched and buffered to a segmented arithmetic unit. This 
segmentation is analogous to an assembly line operation in that some 
elements of the vector will have completed processing at the same time 
that other elements are still undergoing the same processing. The 
effect is that all elementa of the vector appear to have completed 
processing in parallel, hence the term "parallel processor." The STAR 
has been shown to achieve a performance improvement over a CDC 6600 of 
60 to 1 using vector operations. 

The virtual memory system of the STAR-100 gives the illusion that 
physical memory is much larger than it actually is. The portion of the 
user's program not residing in physical memory Is written on a disk in 
the form of 512 or 65,536 64-bit word blocks called pages. When a 
portion of a program residing on the disk is 'needed in main memory, 
the page containing the needed information is moved into main memory and 
the unneeded page is moved out to the disk. Thus, program efficiency is 
improved by structuring code to limit the number of pages moving in and 
out of memory. For a more detailed explanation of paging concepts, 
see reference 4. 

SL/1 (ref. 5) is a programming language developed at NASA/LaRC and 
designed specifically for the STAR-100 computer. The vector type 
instructions facilitate the high performance level of the computer. The 



language also has the capability of using 32-*bit (halfword) or 64-bit 
(fullword) operands. The 32-bit vector will process in half the time of 
a 64-bit vector, but has 32 bits less precision per word. The compiler 
for SL/l executes on a CDC 6400 at NASA/LaRC; thus a program can be compiled 
and edited on a conventional sequential computer before transmission to 
the STAR for execution, resulting in a cost- and time-saving feature. 

In the event of an error, the SL/l language provides convenient diagnostics 
for debugging. This a a time-saving and beneficial feature for the programmer. 
Thus, SL/l was chosen for the program presented in this report because of 
the 32-bit vector capability, the ability to compile on the GDC 6400, the 
error-processing features, and the availability of clear and well -structured 
language concepts. The algorithms used and their translation to SL/l are 
described in the following sections. 

SiJftware Methodology 

A line-by-line radiative transfer computer program was needed to 
efficiently perforin the task of data reduction for infrared pollutant 
sensors being developed by NASA/LaRC. The Simulated Morochromatic 
Atmospheric Radiative Transfer version 2 (SMART2) program was written to 
reduce data-reduction cost by using the vector-processing capabilities 
of the STAR-100 computer and the beneficial features of the SL/l 
programming language. This required a restructuring of the conventional 
algorithms used for line-by-line programs on sequential computers in 
order to conform to vector type operations. The modular structure of the 
program was designed to permit modifications of the computational algorithms 
without affecting the program framework. In addition, thorough documentation 
was necessary for program clarity. 

To minimize computation time without a loss in accuracy, certain 
assumptions were made concerning the atmosphere and radiative transfer 
processes. To compute the transmittance from the surface of the Earth to 
a sensor altitude, h, the modeled atmosphere between these points is 
divided into a number of homogeneous layers, i.e., regions within which 
the temperature, total pressure, and concentrations of the primary and 
interfering molecular absorbers are considered to be uniform. The error 
in this approximation may be made as small as desired by subdividing the 
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atmo$phere into a sufficiently large number of layers. Tlie total radiance 
incident on a sensor at altitude h is given by 

E(h) - / E(u)du (1) 

*'Au 

where Au is the spectral bandpass of interest and E(u) is the total 
monochromatic upweliing radiance. 

For a cloud-free, nonscattering atmosphere under local thermodynamic 
equilibrium, the atmospheric radiative transfer equation for the total 
monochromatic upweliing radiant energy, E(w), seen by a nadir viewing 
type of sensor can be written as 

E(u) ■ eCw) N°(o),Tj) T((d,h) 

nO[w,T(z)] 

0 

+ 7 [1 - e(w)] • (cos 0) Hg(tu) 

• T(u,h) • T(o),h') • f(0) (2) 

where c(u3) is the wavenumber-dependent surface emittance and N^(ui,T ) is 
the Planck blackbody function which is dependent on wavenumber and surface 
temperature T or radiating gas temperature at a particular altitude 
T(z). The monochromatic transmittance of the atmosphere between the 
emitting surface z and the altitude of the sensor h is represented by 
T(o),h), and the monochromatic vertical transmission of the entire modeled 
atmosphere is represented by T(u,h'). The solar zenith angle is 0 and 
the wavenumber-dependent sun irradiance at the top of the atmosphere is 
Hg. The Chapman function f(0) (ref. 1), is equal to sec © for 0° _< 0 
< 60° and is equal to the Chapman polynomial for 0 < 60°. The three terms 
on the right-hand side of equation (2) represent, respectively, the 
surface emission of the Earth, the atmospheric emission, and the solar 
reflected energy (see fig. 1). All of these components must be considered 
in the solar-thermal overlap region at 4.6 wm. These terms are represented 
in SMART2 by TOTSURF, SUMRAD, and RADSUN, respectively. Appendix A 
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Figure 1. Interaction of radiation with the atmosphere . 


contains a dictionary o£ SL/1 variables, and figure 2 is an overview of 
the SMART2 program structure. 

\ closer examination of the form 6f equations (1) and (2) indicates 
two very distinct differences between the methodology of SMART2 and 
conventional line-by-line radiative transfer programs. The first 
difference is the order in which the integrals over Am and over the 
change in altitude, Ah, are performed, By initially Integrating over 
Ah monochromatlcally, the error associated with transmittance averaging 
is eliminated, i.e., the monochromatic transmission at any altitude h 
is given by 

t(u,h) » n t(w,z) C3) 


and not by 

1 ^ 

~(Au,h) • jr T, "rCuiZ) 

z*l 


( 4 ) 


The spectral resolution of the transmittance values, i.e. the number of 
us at which the absorption coefficient is calculated, is defined by a 
constant stepping size, DW. 

The second difference between SMART2 and most other line-by-line 
radiative transfer programs is the specification of the order of limits 
of integration over altitude. The lower limit of integration is the 
radiating source, and the upper limit is the sensor altitude which allows 
for practical computation of the monochromatic and total integrated 
transmission at the top of each atmospheric layer as seen by equation (3) . 
This approach eliminates redundant calculations of atmospheric transmit- 
tance in evaluating signals from aircraft platform sensors at various 
altitudes. 

Theoretically, the total absorption coefficient at any wavenumber 
u consists of contributions from all spectral lines; however, for 





practical computational purposes, only lines within the vicinity of u 
are considered for calculation in SMART2. The contribution to the 
absorption coefficient of lines in the vicinity of u can be divided 
into two parts, direct and wing. Those lines within an interval defined 
by approximately 100 halfwidths of the primary gas to each side of oi 
are considered to be direct contributors to the absorption coefficient, 
while those lines lying outside of this interval result in wing absorption. 
In the SMART2 calculation of direct-line absorption contribution, the 
OMEGA Vector contains center wavenumbers calculated with a constant DW 
stepping size. The spectral lines included for direct- line absorption 
contribution for all elements of the OMEGA vector are determined in 
procedure DETERMN by choosing those lines in the interval OMEGAj l) -CMINV 
to OMEGA|last[ * CMINV. Thus, the vector-processing capability of STAR 
calculates absorption coefficients of all elements of the OMEGA vector 
concurrently with the same spectral lines (WAVEN) . When nev; elements of 
OMEGA vector are calculated, a different WAVEN is computed. The size of 
the OMEGA vector is determined by the user. For atmospheric carbon 
monoxide using the Lorentzian line profile, a SMINV of 5 cm*^ was selected. 
This value is approximately 100 times the line halfwidth of CO, A 250 
element size of the OMEGA vector is sufficiently large to take advantage of 
the vector processing of STAR, but does not affect the direct-line contri- 
bution assumption of a symmetric 5 cm"^ interval about the centerline. 
Presently, the SMART2 program does no', include wing contribution. 

The SMART2 program performs initialization of all variables required 
for the calculation of the monochromatic gas transmissions for center 
wavenumbers in the OMEGA vector (TAU) for each atmospheric layer. The 
algorithm for the computation of TAU is initiated by a call of the READTP 
procedure, which reads all spectral reference information into vectors, 
i.e. line position (OMEGSTR), halfwidth (ALPHA), intensity (SZ) , lower 
energy level (EL), and species identification number (ILINE) from the 
spectral line parameter reference tape for the band interval WI - CMINV 
to WF + CMINV. READTP uses the Q30PNMAP library routine, which performs 
the implicit input by mapping the spectral line file into a COMMON block 
(ABLOCK) , The individual vectors are initialized by equivalencing 
elements of the vectors to elements of the ABLOCK COMMON block. The 
McClatchey spectral line parameter tape (ref. 6) was used to obtain all 
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spectral reference information. The parameters read rhould be under room 
temperature and standard pressure (296 K and 1 atm) conditions or the 
appropriate conversion must be performed prior to reading these parameters 
in procedure READTP (see Appendix B). 

The gaseous transmittance at a particular altitude as a function of 
wavenumber is given by 

TAU(u)) ■ t(w) » enu K. (aj)pc^Jlj (5) 

where Kj(w) is the wavenumber-dependent absorption coefficient of gas 
species i (their sum bej.ng ABSCOP), p (PRES) is the total mean 
pressure of layer k, c^ (GASCONC) is the concentration of absorbing 
gas i, and S. (THICK) is the thickness of the kth layer. 

In general, the absorption coefficient of the nth line of species i 
is described by 




( 6 ) 


where is the layer temperature (TEMP) corrected line intensity and 

3. is the line shape function, i.e., Lorentz, Voigt, or Doppler. The 
Lorentz line shape is given by 


3. (w) 



(7) 


where (ADJALPH) is the temperature- and pressure-dependent line half- 

width for layer k and is given by 


ADJALPH * a. 

in 


“oPe 



( 8 ) 


The calculation of ADJALPH is a vector operation and includes all spectral 
lines in the OMEGSTR vector. The reference line halfwidth, ag (ALPHA), 
is read from the spectral line parameter tape at temperature Tg (REFTEMP) . 
The equivalent pressure, p^, is a function of the ratio of self -broadening 
to the nitrogen -broadening efficiency (BROAD) , the total pressure (PRES) , 
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and the concentration of the absorbing gas (GASCONC) as given by 

■ [GASCONC * (BROAD - 1) + 1] * PRES (9) 

In the case of trace gases in the atmosphere, the equivalent pressure is 
set equal £o the total atmospheric pressure since self -broadening is 
insignificant. 

The line intensity (S) depends upon the temperature through the 
Boltzman distribution factor (ref. 7) and is expressed as 



0‘t 


S » SZ * TCONSQ * EXP | -C8 * (TEMPCON - 1) * EL/REFTEMP| (11) 

where ra is 1.5 for water vapor, methane, and ozone and 1.0 for other 
infrared active molecules, such as carbon monoxide, carbon dioxide, 
nitrous oxide, and C8 is the Boltzman distribution constant; S is 
calculated with a vector instruction in SMART2 and represents the line 
intensities for each spectral line in the OMEGSTR vector. The user has 
the option of calculating ABSCOF as determined by the Lorentz function, 
by the Doppler function, or by an approximation to the Voigt function 
detailed by Pierluissi (ref. 8) in procedure VOIGT. The line shape for 
each layer (IPROF) must be designated in the input and need not be the 
same for each layer. 

Tho total absorption (ABSCOFT) at each wavenumber in the OMEGA 
vector is determined by the summation of the absorption coefficients i‘or 
all lines within the interval. In addition, the user can automatically 
obtain the total absorption coefficient (ABST) that corresponds to a 
maximum of 10 different, primary gas, vertical mixing ratio profiles as 
defined by KMULT. Each value of KMULT is multiplied by the volume 
mixing ratio of the primary gas in every layer, resulting in a bias 
shifting of the input vertical mixing ratio profile. ABST is calculated 
as the sum of two components, ABSCOF(l), which is the absorption coefficient 
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fov the primary gas, and ABSCOFT, which is the absorption coefficient for 
all other interfering gases, The absorption effects of continue , such as 
nitrogen and water vapor, could be easily considered in the form of an 
added third term to ABST after the necessary algorithms for accurately 
calculating the continuum absorption in the spectral region under 
consideration have been developed. 

After completion of the monochromatic calculation of ABST at all 
wavenumbers in the OMEGA vector, the transmittance is calculated as a 
vector instruction shown by equation (5), This procedure is repeated for 
the whole band, resulting in a TAU value for each layer k. The self- 
emission (RADATM) of each layer is then calculated as a vector instruction 
via the wavenumber-dependent energy described by the Planck blackbody 
function. The temperature used in the Planck function is the mean 
temperature (TEMP) of the emitting layer. The emissivity (EMISS) of 
the layer is determined by the calculated transmittance of that layer as 
given by 

EMISS * 1 - TAU (12) 

The total monochromatic transmittance (ATMTAU) at the top of each 
atmospheric layer is then calculated as a vector instruction by 

ATMTAU (uj) = TA,Uj^(ai) (13) 

The next operation performed by SMART2 is the determination of the 
new center wavenumbers to be placed in the OMEGA vector for the purpose 
of repeating the transmittance calculations for the entire spectral band 
being considered, i.e., WI - CMINV to WF + CMINV. The size of the 
wavenumber increment used in absorption coefficient calculations is 
arbitrary and depends upon the spectral resolution required by the user, 
since the absorption coefficient varies rapidly as a function of relative 
position to the spectral line centers. For applications of SMART2 to 
carbon monoxide under tropospheric pressures, a constant stepping size 
(DW) of 0.01 cm"^ was determined to most efficiently describe the wave- 
number-dependent absorption coefficient. The transmittance calculation 
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for a new Oi«IEGA vector is performed and the corresponding results, i,e., 
ATMTAU(w) and RADATM(u), are calculated and summed. This procedure is 
repeated until A'DfTAU and RADATM for the entire ipectral band are 
calculated. 

A program option allows the calculation of the monochromatic emission 
of the Earth (TOTSURF) as a function of one or more surface temperatures 
and emissivlties. This option employs procedure SURFACE in the calculation 
of the wavenumber-dependent Planck distribution of energy as a vector 
operation. Another option to the program, the solar contribution (RADSUN) , 
is calculated as a vector operation in accordance with equation (2) by 
procedure SOLAR, where the Chapman function calculated in procedure 
CHAPMAN is employed in determining the slant path transmittance of incident 
solar radiation. The solar contribution is a function of one or more 
solar zenith angles and surface emissivities. The three terms of equation 
(2), i.e., TOTSURF, SUMRAD, and RADSUN, are then obtained by attenuating 
each component by the appropriate atmospheric transmittance (ATMTAU) 
value for an altitude corresponding to the top of each atmospheric layer. 
For convenience, the trapezoidal rule is used for integration of all 
variables over the spectral band being examined. Computed results are 
then printed out. A complete program listing and brief description of 
each of the program sections are contained in Appendix C. 

Operating Instructions 

The input from cards contains the atmospheric, surface, and solar 
parameters used by SMART2. The physical setup for the PRES, TEMP, THICK, 
and IPROF data is arranged such that the atmospheric layers are read in 
ascending order of altitude (descending order of pressure) . The GASCONC 
data follow in the same order with the concentrations of the gases in 
each layer placed in the 10-column field corresponding to the identi- 
fication number of the gas, i.e., water vapor (ID » 1) concentration in 
columns 1 to 10 and carbon dioxide (ID = 2) in columns II to 20. 

To avoid unnecessary calculations, one should exercise care in 
choosing the values of DW and CMINV. For our sample case using 
carbon monoxide as the primary gas, 2080 to 2170 cm ^ as the band width, 

DW of 0.01 cm'^, and CMINV of 5 cm~^, 90,000 points of integration were 
considered with 1,557 spectral I’.r.es included from the tape. If DW 



were chosen smaller and/or CMINV were chosen larger, these totals would 
be greater, causing execution time to increase, 

The size of the OMEGA vector is determined by the LITERALLY declaration 
of NINC in MODULE SMART2. Since the STAR is a vector-processing machine, 
use of vector lengths greater than 100 is recommended. The vector size 
should be sufficiently largo to take advantage of the machine's hardware, 
but should not affect the 5 cm"'^ assumption on the direct conti’ibution 
interval about the centerline. The optimum size of OMEGA was determined 
to be 250. 

Caution should be exercised concerning the units of the line para- 
meters read from the spectral line tape. These units must be as specified 
in the dictionary of SL/1 variables, and the reference temperature of 
the spectral lines must be 296.0 K or the appropriate change to REFTEMP 
in the program must be performed. A sample case input is listed in 
Appendix D. 

Presently, one may choose either Lorentzian, Doppler, or Pierluissi's 
approximation to the Voigt profile in the calculation of spectral line 
shape. The Voigt profile requires a significant increase in the calcu- 
lation time. 

The lOPT parameter allows the user to choose which components of 
the equation of radiative transfer [eq, (2)] will be calculated. An 
lOPT of 1 will calculate only the atmospheric portion for 1 atmospheric 
profile, up to 10 concentrations of the primary gas using the XMULT 
vector, and up to 9 different interferents. An lOPT of 2 will calculate 
the atmospheric contribution and the Earth surface omission for one or 
more surface temperatures (or surface emissivities) . An lOPT of 3 will 
calculate the atmospheric and surface contributions in addition to the 
solar reflected energy for one or more solar zenith angles and surface 
emissivities. 

The program is currently set up to perform carbon monoxide calcu- 
lations in the 4.6-pm spectral band. Any other primary pollutant gas 
may be considered by setting the values of IPOLLUT and IDENT(l) to 
the identification number of that species and by making necessary 
adjustments to spectral band widths. The gas-broadening coefficient, 

BROAD, must be changed to correspond to the gas being considered. 
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The trapezoidal rule is used as the integration approximation. Should 
the user wish to incorporate a different integration process « the trapezoidal 
rule may be easily replaced, 

The program output consists of the band-integrated average transmission 
and total upwelling radiance as a function of altitude and primary gas 
concentration. In addition, the band-integrated total surface radiance 

t 

is listed as a function of altitude, surface temperature, surface emis- 
sivity, and primary gas concentration. The band -integrated total solar 
radiance is printed as a function of altitude, solar zenith angle, surface 
emissivity, and primary gas concentration. A sample case output is listed 
in Appendix E. 

The input listed in Appendix D and output in Appendix II correspond 
to a problem of determining the atmospheric, surface, and solar radiance 
calculated at the top of each of the 15 modeled atmospheric layers. The 
model chosen is a midlatitude summer atmosphere with a corresponding 
water vapor profile (ref. 9). The sun zenith angles are 45® and 75®, The 
pollutant is carbon monoxide, calculated at 10 different concentrations in 
each layer, with interferent constituent concentrations of water vapor 
and carbon dioxide held constant in each layer, The spectral band is 
4.6 ym. 


Comparisons and Conclusions 

Any computer program is only as accurate as the theoretical model on 
which it is based and the numerical algorithms coded into program 
instructions. As outlined previously under "Software Methodology," the 
accuracy of this model is dictated by the assumptions made in our 
application to a specific problem. 

To verify the results of the SMART2 calculated absorption coefficients, 
a comparison was made to an existing program written by Casas and Campbell 
(ref. 2). This program has been tested and verified with analytical 
calculations of absorption and laboratory spectra. The comparison was 
made on the basis of computation time and accuracy, For a 15- layer 
atmospheric model including carbon monoxide, water vapor, and carbon 
dioxide, for a 90 cm’^ band interval, at a constant stepping size of 0.01 
cm"^, using 2 surface temperatures, 2 surface emissivities, and 2 solar 
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zenith angles, assuming a Lorentzian profile in each layer, the total 
processing time was 2898,79 sec on a Cyber 173. On the STAR-100, the 
processing time was based upon the size of the OMEGA vector. The optimum 
vector size was determined by comparing four test cases using vector 
sizes of 250, 500, 1000, and 2250 for the atmospheric contribution only. 

The processing times were (a) 250 vector size - 16.58 sec, (b) 500 

vector size - 16,08 sec, (c) 1000 vector size - 14.61 sec, and (d) 2250 

vector size - 13.89 sec. A comparison was made between the band-integrated 
average transmission of the Casas and Campbell program and the SMART2 
program with the four vector lengths. The results are shown in tables 
1 and 2. Thus, a vector size of 250 is the most accurate. Using this 
vector size, the same test case run on the Cyber 173 was run on the STAR-100, 
resulting in an execution time of 29.2 sec. This is a factor of 99 improve- 
ment over the Cyber 173. 

The processing times for 11 test cases using the same 15 layer 
atmosphere run on the STAR computer are listed below. Although a vector 
length of 250 has been shown to be the most accurate, a vector length of 
1000 was used for these test cases to save time and cost. 

SL/1 has an optimizing parameter (OPT) associated with the compiler 
that is used to optimize code. The higher the optimization parameter, 
the more efficient the code. The test cases run for vector lengths of 
100 are listed in table 2. Thus, the Voigt profile takes significantly 
more execuvion time than the Lorentz. The surface contribution adds an 
insignificant amount of time, while the solar calculation adds about 2 
sec per case. 

A line-by-line radiative transfer program for the STAR-100 computer 
has been developed and described. Any of the specific task procedures may 
be easily substituted by algorithms that are more suitable and convenient 
to the user. 

The program currently performs carbon monoxide calculations in the 
fundamental 4,6-um spectral band. By making necessary adjustments. 
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Table 1. 


Percentage difference between band-integrated average transmission 
of ^1ART2 program and Casas and Campbell program. 


Vector 


Concentration of CO fooml 

Size 

Layer 

0.0 

0.2 

l.O 

250 

1 

0.081 

0.088 

0.099 


5 

0.206 

0.229 

0.300 


10 

0.221 

0.255 

0.364 


15 

0.221 

0,258 

0.384 

500 

1 

0.123 

0.131 

0.151 


5 

0.322 

0.352 

0.4S3 


10 

0.348 

0.394 

0.555 


15 

0.349 

0.399 

0.587 

1000 

1 

0.185 

0.196 

0.224 


5 

0.492 

0.537 

0,683 


10 

0.533 

0.601 

0.83S 


15 

0.534 

0.609 

0.878 

2250 

1 

0,231 

0.245 

0.286 


5 

0.621 

0.679 

0.880 


10 

0.673 

0.763 

1.086 


15 

0.674 

0.776 

1.144 
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Table 2. Execution times for I I test cases run on the STAR- 100, 


Test Case 


1. Atmospheric contribution only, all layers using Lorentz 
profile 

a. OPT » 1 

b. OPT ■ 2 

c. OPT » 3 

2. Atmospheric contribution only, Doppler profile used in 
layers 14 and 15 

a. OPT - 1 

3. Atmospheric contribution only, all layers using Voigt 
profile 

a. OPT - 1 

b. OPT « 2 

4. Atmospheric contribution, surface contribution with 2 surface 
temperatures and 2 emissivities, all layers using Lorentz 
profile 

a. OPT « 1 

b. OPT « 2 

5. Atmospheric contribution, surface contribution with 2 surface 
temperatures and 2 emissivities, solar contribution with 2 
solar zenith angles, all layers using Lorentz profile 

a. OPT » 1 

b. OPT * 2 

6. Atmospheric contribution, surface contribution with 2 surface 
temperatures and 2 emissivities, solar contribution with 2 
solar zenith angles, all layers using Voigt profile 

a, OPT » 2 


Time (sec) 


14.61 

12.48 

12.46 


20.62 


119.10 

118.89 


14.74 

12,79 


24.28 

21.80 


131.75 
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other gas species transmittance calculations may be made, Future additions 
to the program will include a sub-Lorontzinn wing absorption algorithm, 
the water vapor and nitrogen continuum absorption algorithms, and a routine 
to plot transmittance as a function of wavenumber. 

The storage required for the program is not applicable to the STAR-100 
computer, For our sample case, 15 layers required 29,2 sec of execution 
time. By utilizing the vector capability of the STAR computer, program 
efficiency incrf'ased, and a cost savings factor of seven was realized. 

ST. PETERSBURG FLIGHT TEST 

During the period of August 13 to August 20, 1979, a series of flight 
tests of the MAPS (Measurement of Air Pollution from Satellite) brass- 
board sensor system test evaluations was performed in the St. Petersburg, 
Florida area. The purpose of this series of flight tests was to provide 
a data base for the evaluation of the response of the carbon monoxide 
channels of the sensor system to (1) hot and humid subtropical atmospheric 
conditions, (2) partial cloud cover in the IFOV of the sensor, and (3) land/ 
water interface changes. 

The urban area of St, Potersburg-Clearwater provided a relatively 
isolated site from industrial carbon monoxide emission. The largest 
production of CC in this area results from automobile traffic. An estimate 
of the emissions was obtained by correlating traffic counter data with EPA 
automobile emission rates. 

The remote sensor system and supporting equipment were flown onboard 
a Cessna 402B twin-engine aircraft modified with a nadir viewing port. The 
aircraft was equipped with dual ?60 channel navcoms, DME, ADF, and a Loran 
navigation system. The aircraft experimental hardware included (1) remote 
carbon monoxide sensor, (2) radiation thermometer (11.6 urn), (3) hygro- 
meter system, (4) total air temperature probe, (5) altitude pressure 
transducer, (6) LORAN route verification system, (7) data acquisition 
system, (8) 35-mm aerial camera, and (9) LORAN navigation system. 
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The test was comprised of one primary test area and three secondary 
test areas* The primary area was the St. Petersburg-Clearwater peninsula 
as shown in figure 3. Test flights at an altitude of 6.1 km (20,000 ft) 
and 2,4 km (8,000 ft) were flown during the test period. An extensive 
report of the results of these tests is forthcoming. In general, the flight 
test was a success and an excellent data base was obtained. A work task 
involving a MAPS commitment to be flown on the Space Shuttle delayed the 
reduction and analysis of the data from the St. Petersburg flight; 
however, a complete analysis of it will be performed prior to December 
1981. 


SUNWER MONEX EXPERIMENT 

The Summer Monsoon Experiment (MONEX) was conducted over the Indian 
Ocean from May to August 1979 in conjunction with the second Special 
Observing Period of the Global Weather Experiment. The MONEX observations 
provided an unprecedented data set for a basic study of the monsoon 
phenomenon. The Asian monsoons, due to their dimension and intensity, 
represent the largest disturbance of general circulation, creating strong 
cross-equatorial flows and hemispheric interactions, During this test 
period the MAPS remote sensor was operational in over 20 test flights. The 
purpose of this research was to examine the Interhemispheric exchange of 
CO during the monsoon. 

ODU personnel performed research in two primary areas supporting 
NASA/LaRr. work. The first area of research support was in the operation 
of a gas chromatographic system which obtained CO and CHi, concentration 
profiles from the CV990 aircraft during the test period. A very large 
data base was obtained. As collaborative scientists, ODU personnel 
performed research tasks in the planning of the tost, new software 
development and testing of a new inflight data reduction algorithm. 

Specific tasks performed involved extensive analysis of meteorological 
data and radiometer data obtained from the MAPS remote sensor system. 

High altitude (12-km) remote sensor data was obtained for the first 
time. An extensive data set was acquired during the Saudi Arabia 
Experiment and Arabian Sea Experiment phase of the Summer MONEX program. 




Figure 3. Pri»ary test ar; 





A comprehensive analysis of the entire data set and publication of the 
results will be completed by October 1981. 

OSTA-1 PAYLOAD 

Investigative studies were implemented for the development of a MAPS 

f 

©xpc'iment for the OSTA-1 payload for §pace Shuttle, Comprehensive 
parametric studies were utilised in the preliminary design and development 
of this space platform version of the MAPS remote sensor system, 

These theoretical studies were required to determine operational 
limitations for a MAPS global remote sensor CO measurement experiment, The 
implementation of the MAPS experiment to a space platform involved many 
multifaceted fundamental research studies and involved many man hours. 
Specific areas of work were optimisation of instrument gas cell pressures 
and determination of (1) Instrument signal sensitivity to atmospheric 
temperature profiles changes, (2) instrument signal sensitivity to 
surface temperature characteristics, (3) gas filter correlation instrument 
signal function, and (4) calibration requirements. 
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APPENDIX A 


DICTIONARY OF SL/1 VARIABLES 


MDOULE SMART2 

A 

AA 

ABSCOF(fiS,NL)lNINCl 

abscoft|ninc| 

abst|ninc| 

adjalph|nlin1 

alpi-ia|nlin| 

alphd|nlin| 

ATMTAU(NX)|NINN| 

bbInlinI 

bd|ninc| 

be|ninc| 

bf|ninc] 

bradx 

brdfac 

broad 

C7 


number of integration points 

length of wavenumber interval, cm’^ 

monochromatic absorption coefficients 

monochromatic total absorption coefficient 
for all interfering gaseous constituents 

monochromatic total absorption coefficient 
for all gaseous constituents 

adjusted halfwidths for all spectral lines 
considered 

halfwidth values at half maximum for all 
spectral lines considered 

halfwidth values at half maximum for all 
Doppler spectral lines considered 

calculated monochromatic transmission of the 
atmosphere at the top of the atmospheric 
layer 

bit vector used to find the nonlinear 
molecular species for rotational partition 
function calculation 

bit vector used to determine if REGIONl 
of Voigt profile should be used 

bit vector used to determine if REGI0N2 
of Voigt profile should be used 

bit vector used to determine if REGIONS 
of Voigt profile should be used 

primary gas pressure broadening coefficient 

primary gas pressure broadening factor 

short form of BRADX 

first radiation constant (1.1908E-12 W cm"^ 
sr"^cm"^) 
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C8 

CHAP N2 
CM2 


Boltzmann's constant (1.439 cm K) 

exponential component of the Chapman transmis- 
sion function of the atmosphere for all solar 
zenith angles 

2 * CMINV 


CMINV 

C0EFl9| 

COMPABSjNINCl 


miniimim distance above and below OMEGA 
considered as centerline absorption in the 
transmittance calculation 

coefficients of eighth degree polynominal 
used in Chapman function 

line profile component of absorption 
coefficient 


conc|ns 


concentration of gases used for transmittance 
calculations 


D1 


0 . 


D2 

DIFSUN[411 

DIFW|4lj 

dunwy(nl)1nlin1 

DUMMYI(NL)|NLIN1 

DW 

DWl 

le|nlin| 

emiss|nem| 

emisxInemI 

FIRST 

GASCONC (NL)|i6l 


1 . 

vector of differences between succeeding 
elements of SUNFLUX vector 

vector of differences between succeeding 
elements of SUNW vector 

dummy vector tjsed for intermediate calculations 
dummy vector used for intermediate calculations 
integrating wavenumber increment, cm”^ 

32-bit form of DW 

energy of lower state of transition of 
OMEGSTR from spectral line parameter tape, 
cm"^ 

32-bit form of EMISX 

thermal emissivities of the surface of the 
Earth 

the position of the first spectral line to 
be included for centerline absorption 
relative to the OMEGSTR vector 

32-bit form of GASCONX 
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GASCONX(NL) 16 
I 

IDENTl 10 i 

identn|nit| 

iline|nlin| 

lOPT 


IPOLLUT 

iprof|nl| 


lUPPER 

J 

K 

L 

LOWER 


concentration of gases used for transmittance 
calculations, ppm 

do loop inde;c 

vector of identification numbers corresponding 
to gas species identification numbers from 
spectral line parameter tape 

vector of gas species identification numbers 
for spectral lines to be included in center- 
line absorption 

vector of gas species identification numbers 
of OMEGSTR from spectral line parameter tape 

option parameter 

1 - atmospheric contribution only 

2 - atmospheric and surfat?'^ contributions 

3 - atmospheric, surface, and solar 

contributions 

primary gas species identification number 

line profile option parameter 

1 - Lorentz 

2 - Voigt 

3 - Doppler 

number of integration points calculated 
do loop index 
do loop index 
do loop index 

lowest wavenumber in OMEGSTR vector 


M do loop index 

MCCLAT] 12800| vector of spectral line information placed 

on Q30PNMAP file 

M0LWT|NLIN| molecular weights of gas species corresponding 

to OMEGSTR vector 

MWT|12| molecular weights of gas species corresponding 

to gas species identification number 

NEM number of surface emissivities used for 

dimension purposes 
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NEMIS 

number of actual surface emisslvitles 
considered 

NINC 

number of OMEGAs considered for transmittance 
calculation at a time 

NIT 

maximum number of wavenumbers considered 
for centerline absorption calculation used 
for dimension 

NL 

number of atmospheric layers used for 
dimension purposes 

NL\Y 

number of actual atmospheric layers considered 

NLIN 

number of spectral lines from tape used for 
dimension purposes 

NaMEG 

number of actual spectral lines from tape 
considered 

NS 

number of gas species used for dimension 
purposes 

NSG 

number of surface temperatures used for 
dimension purposes 

NSPEC 

number of actual gas species considered 

NSURF 

number of actual surface temperatures 
considered 

NUMN 

number of actual wavenumbers considered for 
centerline absorption calculation 

NX 

number of primary gas concentration multipliers 
used for dimension purposes 

NXM 

number of actual primary gas concentration 
multipliers considered 

NZ 

number of solar zenith angles used for dimension 
purposes 

NZEN 

number of actual solar zenith angles considered 

omega] NINC 1 

vector of center wavenumbers of the sub- 
intervals being considered for transmittance 
calculations 

OMEGSTRjNLINj 

vector of wavenumbers read from spectral line 
parameter tape 
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OPATHI NI„1 

optical path of atmospheric la>^ers, atm-cm 

pInincI 

vector of Ioj-oj ! 

1 o' 

planck|ninc| 

thermal emission of the surface 

pres|nl| 

32 -bit form of PREX 

prex|nl| 

pressure of atmospheric la>^ers, atm 

PR0FAC1(NL)|NLIN| 

multiplicative variable for line profile 

PROFAC2iNINC| 

multiplicative variable for line profile 

C|w-Wq|/o3 

RADATM(NX)|NINC| 

the radiation upwelling from the atmosphere 
as a result of atmospheric gaseous molecular 
emission 

radsun|nx| 

SUMSUN integrated 

REFTEMP 

reference temperature of spectral line 
parameters from tape 

rootemp|nl| 

square root of layer temperature 

S(NL) |NLIN| 

adjusted line strength of spectral lines 
from tape 

sumrad(nl3|nx| 

RADATM integrated 

SUMSUN(NL,NZ,NEM) |NX| 

radiation upwelling from the atmosphere as 
a result of incident solar energy 

SUMSURF (NL , NSF , NEM) | NX | 

radiation upwelling from the atmosphere as 
a result of Earth surface emission 

SUMTAU(NL)|NX| 

ATMTAU integrated 

sun|ninc| 

interpolated solar flux values for OMEGA 
considered 

SUNC0M(NEM)|NZ| 

wavenumber -independent portion of solar 
contribution 

SUNFLUX|4l| 

wavenumber-dependent solar energy incident 
at the top of the atmosphere 

SUNW|4l| 

wavenumbers corresponding to SUNFLUX 

surfcom|ninc| 

component of surface energy contribution 

sz|nlin] 

line strength values from spectral line 
parameter tape, atm ^cm ^ 
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f 


tau|ninc| 

transmission of an atmospheric layer 

TAUSLAT(NX,N2)|NINCl 

total atmospheric slant path coefficient 

TCONSQ(NL)[NLIN| 

temperature dependence of the rotational 
partition function 

TEMPI NLI 

32-bit form of TEMX 

TEMPCON(NL)iNLIN| 

temperature adjustment for line halfwidth 

temx|nl1 

temperature of atmospheric layer, Kelvin 

THETA 1nz| 

solar zenith angles 

THICK 1NL| 

32-bit form of THICX 

thicxInlI 

thickness of atmospheric layer, cm 

tinv|nl| 

inverse of layer temperature 

totsunInincI 

radiation upwelling from the atmosphere as 
a result of incident solar energy 

TOTSURFInxI 

SUMSURF integrated 

TRAN (NX, NL) In INC 1 

vector used to save atmospheric transmission 
values 

TSURF 1 NSF 1 

32 -bit form of TSURX 

tsurxInsfI 

surface temperatures to be considered 

UPPER 

highest wavenumber in OMEGSTR vector 

Ifl 

the first wavenumber of the OMEGA vector 

W3lNINC| 

OMEGA cubed 

wavenInitI 

vector of wavenumbers to be included in 
centerline absorption 

WF 

final wavenumber used for transmittance 
calculation 

WI 

initial wavenumber used for transmittance 
calculation 

XjNINCl 

Voigt absorption coefficient parameter 

XMULTjNXj 

32-bit form of XMULX 

XMULXlNXj 

primary gas concentration multipliers 
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y NINC 


Voigt absorption coefficient parameter 


PROCEDURE READTP 
lERR 

PROCEDURE DETERMN 
BA 

IJ 

OMEGl 

C»fEG2 

PROCEDURE REGIONl 
AN|301 

Cl, SI, S3, T, T3, XI, XSER, 

Yl, yser|ninc| 

LL 

N, NCOUNT 

PISQ, X2, XN, YNEW, Y2, YN, 
YNEW 

PROCEDURE REGION2 


Al, 

A2, A3, 

A4, 

AS, 

A6 

Cl, 

F, G, H 

, R. 

SI 

S3, 

T3, 

XI, X2, 

Yl. 

Y2| 

NINC 


PROCEDURE REGION 3 

Bl, B2, B3, B4 

Cl, F, G, R, SI, S3, T, T3, 
XI, Yl|NINCl 


error parameter for Q30PNMAP 

bit vector used to determine the lines to be 
included in the interval about the center 
wavenumber for transmittance calculation 

do loop index 

first wavenumber in the interval about the 
center wavenumber 

final wavenumber in the interval about the 
center wavenumber 

Voigt profile coefficients 

Voigt profile parameters 

do loop index 

Voigt profile parameters 

Voigt profile parameters 

Voigt profile parameters 
Voigt profile parameters 

Voigt profile parameters 
Voigt profile parameters 
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APPENDIX B 


UNIT CONVERSION 


In the past, various units have been used to define the line intensities; 
therefore, the following conversion factors may be hejnful. 

At standard temperature and pressure condition, STP, 

1 Ccm-atm)g^p ■ 2.69 x 10^ ^ molecules/cm^; 
however, at some temperature T', 


1 (cm-atm).p, 


X 2.69 X I0l9 molecules/cm^ 
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APPENDIX C 


PROGRAM LISTING 
Module SMART2 

SMART2, the program executive^ defines the global declarations. The 
LITERALLY declaration associates an Identifier with a sequence of characters. 
The INITIAL declaration initializes variables prior to execution. 



a 

II 


• 5 


m 

at 
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Procedure VARIABLE 


VARIABLE is a procedure subordinate to SMART2, VARIABLE initializes 
variables to appropriate values. 
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Procedure TRANSMT 


TRANSMT is a procedure called by VARIABLE. TRANSMTT performs the 
following*. 

(1) Reads the atmospheric profile information by calling READTP; 

(2) Calculates the temperature- and pressure-dependent line half- 

I 

width and the line intensity for all spectral lines j 

C33 Determines the wavenumbers at which the absorption coefficient 
is calculated using a constant stepping size; 

(4) Chooses the spectral lines included in the direct line 
absorption contribution by calling DETERMNj 

(5) Calculates absorption coefficients using either the Lorentz, 
Doppler, or Pierluissi's approximation to the Voigt profile; 

C6) Calculates atmospheric transmittance and radiance as a function 
of altitude and primary gas concentration; 

(7) Includes surface monochromatic emission as a function of one or 
more surface temperatur*es and emissivities; and 

(8) Calculates solar contribution as a function of one or more 
solar zenith angles and surface emissivities. 

Functions (7) and (8) are program options. 
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Procedure READTP 


READTP is called by TRANSMT. READTP reads the spectral line parameter 
file for all lines from WI - CMINV to WF + CMINV. The Q30PNMAP library 
routine performs implicit input by mapping the spectral line file into the 
ABLOCK COMMON block • The individual vector elements are initialized by 
equivalencing to elements of the ABLOCK COMMON block. 
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Procedure DETER^tN 


OETEWfN is called by TRANSMT. DETERMN determines the spectral lines 
included about the center wavenumber vector (OMEGA), The lines included 
fall into the interval OMEGA(l) - CMINV to OMEGA(last) + CMINV. 
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Procedure VOIGT 


VOIGT is called by TRANSMT. VOIGT determines which region of Pierluissi s 
approximation for the Voigt profile should be applied. 
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Procedure RBGIONl 


REGIONl is called by VOIGT. REGIONl calculates the Voigt absorption 
coefficient for X and Y in region 1 (see ref. 8) of Pierluissi's method. 
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Procedure REGION2 


REGI0N2 is called by VOIGT. REGION2 calculates the Voigt absorption 
coefficient for X and Y in region 2 (see ref. 8) of Pierluissi's method. 
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Procedure REGIONS 


REGIONS is called by VOIGT, REGIONS calculates the Voigt absorption 
coefficient for X and Y in region S (see ref. 8) of Pierluissi's method. 
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Procedure SURFACE 


SURFACE is called by TRANSMT. SURFACE calculates unattenuatod surface 
emission using the Planck blackbody function. 
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Procedure CHAPMAN 


CHAPMAN 

transmission 


is called by TRANS^r^. CHAPMAN determines the slant 
using the Chapman polynomial (ref, 1). 
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Procedure SOLAR 


SOLAR i» called by TRANSMT. SOLAR calculates unattenuated solar 
radiation reflected by the surface of the Earth. 
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APPENDIX D 


INPUT PARAMETERS FOR A SAMPLE 
TEST CASE LISTED AS CARD IMAGES 
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APPENDIX E 


OUTPUT LISTING FOR A SAMPLE TEST CASE 
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